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We start with quasi-exactly solvable (QES) Hermitian (and henee real) as well 
as complex PT-invariant, double sinh-Gordon potential and show that even after 
adding perturbation terms, the resulting potentials, in both cases, are still QES 
potentials. Purther, by using anti-isospectral transformations, we obtain Hermitian 
as well as PT-invariant complex QES periodic potentials. We study in detall the 
various properties of the corresponding Bender-Dunne polynomials. 



INTRODUCTION 



In recent years, various features of the complex PT-symmetric Hamiltonians have been 
explored in the literature,see for example Refs. [l-4|. In particular, it has been shown 
that so long as the PT-symmetry is not spontaneously broken, then the energy eigenvalues 
of the Schrodinger equation are real. Further, several PT-symmetric complex QES P, 
potentials have been discovered ¡8l4l0l|. The purpose of this paper is to point out that even 
after suitably perturbing either a Hermitian or a complex PT-invariant QES potential, 
one can still obtain another Hermitian or complex, PT-invariant, QES potential. 



II. PT SYMMETRIC NON-HERMITIAN QES SYSTEM 

A. The Model 

We start from the well known non-Hermitian but PT-invariant potential — [a cosh(2x) — 
íMY , (^ = 1 = 2m) ¡8| and show that even the perturbed complex PT-invariant Hamil- 



tonian 



H = p^-[a cosh(2x) - ^Mf + - (2.1) 

smh [X) cosh [x) 
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where a, /, M are real and — 1 < / < O is a QES potential so long as M — 21 — 1 or M + 2/ + 1 
is a positive even integer. Note, the restriction — 1 < / < O has been imposed so that the 
potential is not too singular at a; = O and there is a communication from the left side to 
the right side. 

Let US first show that the above non-Hermitian Hamiltonian (12. ip is a PT-symmetric 
one. We note that in this case the parity transformation is defined which is 

reflection of coordínate about the point x = ij. Under the time reversal transformation 
t — y —t and further, one replaces i — > —i. One can check easily that, under PT 
cosh(2x) — > — cosh(2a;) ,cosh^(x) — > — sinh^(a;) and sinh^(x) — > — cosh^(x). 

We substitute 

^{x) = e'^^°'^^^''U{x) (2.2) 
in the Schroodinger equation Hip = Eip with the H given by Eq. (12.11) and obtain 

0"(x) + 2msinh(2a;)0'(x) + [{E - + a^) 

-2z(M - l)acosh(2x) - -^i^Ut^ + ^''^ X}\ ]<P{x) = O 

smh [xj cosh [xj 

On further substituting 

(f) = [cosh(a;)]"[sinh(x)]^?7(a;) , (2.3) 
we find that the system admits non-singular QES solutions provided 

a{a + 1) = /(/ + l)and (3{(3 + 1) = /(/ + 1) (2.4) 
is satisfied. The condition in Eq. (12. 4p implies either of the following four conditions i.e. 
(i) a = 13 = 1 + 1 (ii) a = ¡3 = -l (iii) a = = l + 1 (iv) a = l + I,l3 = -l. 
We discuss all the cases separately. 
Case (i): a = /3 = 1 + 1, 

In this case, 7]{x) as given by Eq. (12. 3p can be shown to satisfy 

r]"{x) + 2[{l + 1) coth(a;) + (/ + !) tanh(a;) + ia smh{2x)]r]' (x) 



+ 



e — iza cosh^(xl 



r]{x) = 0, (2.5) 



where 



e = E - + + 4(1 + 1^ - 2ia{2l - Ai + 3) 

z = 4(M-2/-3). (2.6) 
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This system has p QES solutions in case M — 2/ — 1 = 2p with p = 1, 2, 3, ■ ■ ■. These 
Solutions are of the form 

oo 

r]{x) = Y.a,,{i)^[cosh2{x)r, (2.7) 

n=0 

Two of the low lying solutions are 

M = 2/ + 3 , r] = constant , E = -a^ + 4/ + 5 , (2.8) 

, , , A E + a^ -121-21 

M = 21 + 5, 71 = A cosh 2x + iB , — = , 

B 4a 

E = 8l + 15-a^± ^{21 + 3)2 - Aa^ . (2.9) 

Note that since a = ¡3 the PT symmetry is unbroken and henee, as expected, the energy 
eigenvalues are real. 

Case (ii): a = (3 = —l. 

The solutions in this case can be obtained from the solutions of case (i) by everywhere 
changing / — > — / — 1. Thus in this case, one has p QES solutions in case M + 21 + 1 = 2p 
with p = 1, 2, 3, ■ • •. Note that since a = P, even in this case, the PT symmetry is 
unbroken and the energy eigenvalues are real. 

Case (iii): a = —l , /3 = l + 1. 

In this case we have p QES solutions in case M = 2p, p = 1,2, ■ ■ ■ . Two of the low 
lying QES solutions are 

M = 2, E = 3-a^ -2ai{l + 2l), (2.10) 



M = 4, E = ll-a^-2ai{l + 2l)±^{l-a^)-ia{l + 2l). (2.11) 

Note that since in this case a ^ ¡3, henee the PT symmetry is broken spontaneously and 
eigenvalues are no more real. In the special case of a = /3 = |, the PT symmetry is 
restored and E becomes real. 

Case (iv): In this case, the solutions can simply be obtained from the solutions of the 
case (iii) by changing / — > —l — l. 

Ordinarily, the boundary conditions that give the quantized energy levéis are il){x) — > O 
as |a;| — i- oo on the real axis. However, in the present case, we have to continué the 
eigenvalue problem into the complex —x plane^]. On putting, x = u + iv where u and 
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V are real, it is easy to see that for u > O the boundary condition is satisfied so long as 
— TT < f < — I (mod tt) while for n < O it is satisfied if — ^ < i; < O (mod tt) 

It is worth mentioning that, the non-Hermitian PT symmetric Hamihonian in Eq. (12. ip 
after a suitable change of variable can be expressed in terms of the SL{2,R) generators ( 
at most quadratic). To show this we substitute t = cosh(2a;) in the Eq. (12.51) to obtain 
Hgr] = Ef] where, 

Hg = -A{f - 1)^ - [(8/ + 12)í + Aatif - 1)]^ - [-M^ + 

+ 4{l + 1)2 - 2ai{2l - M + 3) - ^(í + 1)](2.12) 

This gauged Hamiltonian then can be expressed in terms of the generators of the SL{2, R) 
by 

Hg = -4[(J2 - J!) + m(J+ - J_) + {n + 2l + 2) Jo] - [-M^ + 

+ 4(/ + l)2 + 4n(/ + l)+n2] (2.13) 

while the generators are given by J_ = ^, -Jo = ~ ^ and 7+ = t^-^ — nt. Gauged 
iamiltonian in terms of SL{2,R) for a more general system has been discussed in Ref. 

3- 



B. Bender-Dunne (BD) Polynomials 

Case (i): a = l + 1 = ¡3. We make a change of variable in Eq. (l2.5p . cosh2(x) = t , 
yielding 

t{t - l)V'(í) - [(/ + 3/2) - (2/ + 3 - 2ia)t - 2mí>'(t) + 

+ {l/A)[e-iazt]ri{t) = (2.14) 

On further substituting 

V{t) = t ^r(^ÍfLfo^ ^ (2-15) 
^on\r{n + 1 + 3/2) 

yields the three-term recursion relation satisfied by the poljTiomials, -Pn(e) 

Pn+i{e) - [^ + n{n + 2l + 2-2ia)]Pnie) 

+ ia{n + l + l/2)n[M - 21 - 2n - l]P„_i(e) = O . (2.16) 
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First few Polynomials are: 

m = ¡ 



e 3 
P2{e) = — + -{21 + 3 - 2ia) + ia{l + -){2l + 3 - M) (2.17) 
lo 4 2 



For the case, M = 2Z + 3, -Pi(e) is the critical Polynomial ll|. On demanding Pi(e) = O 
correctly yields the QES energy eigenvalue as given by Eq. (12. 8p . On the other hand, for 
M = 2¿ + 5, -P2(e) is the critical Polynomial, and demanding P2(e) = O correctly yields 
the QES energy as given by Eq. (12. 9p . 



Following Bender and Dunne 
polynomial. We find 



it is easy to compute the norm (7„) of the n-th 



7„ = {AaiYn\ [| (A; + / + 1/2) (M - 2A; - 2/ - 1) (2. 18) 



fc=i 



Weight factors (wi and uj2 ) [1 1] , for M = 2/ + 5 are 

2/ + 3 - 2az 1 
= , = + - 

^{21 + 3)2 - 4a2 2 

2[+3-2a^ 1 

W2 = , = + o (2-19) 

^(2/ + 3)2 - 4a2 2 

Moments of weight function is defined by 

^ir^ = j dEu{E)E'^ (2.20) 

It is easily shown that the n-th moment, for large n, is proportional to (M + a^)^. 

One can similarly study the properties of the Bender-Dunne polynomials in the other 
three cases. 



III. PERIODIC PT INVARIANT QES SYSTEM 



As has been shown in 
the potential v{x) - 



13|, if under the anti-isospectral transformation x — )■ ix = ?/, 



f (y) and if the potential v{x) has m QES levéis with energy 
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eigenvalue and eigenfunctions E¡¡. (fc = 0, 1, 2 ■ ■ ■ m — 1) and ipkix) respectively then the 
energy eigenvalues of v{y) are given by 

Ek = -Em-l^k, i'kiy) = 1pm-l-k{Íx) (3.1) 

Under this anti-isospectral transformation,x — > ix = 6, it is easily seen that the Hamil- 
tonian ( I2.ip goes over to 

H = p' + [a cosm -^Mr + + ^-^^ , (3.2) 

with — 1 < / < 0. In this case U{x) = U{x + ir). This complex Hamiltonian is invariant 
under combined Parity [ x — > f ~ ^] Time reversal [ t — )■ —t Szi — )■ —i]. Explicit 
QES Solutions can be obtained from the solutions of hyperbolic case discussed in detall 
in section 2.1 by using the anti-isospectral transformation as given in Eq. (13. ip . One can 
easily construct the BD polynomials and their properties for this case by foUowing the 
methods outlined in section 2.2. 



IV. REAL QES SYSTEMS 

Before completing this paper, it may be worthwhile to point out that, starting from 
the Hermitian, QES, DSHG potential V{x) = [acosh(2x) — M]^, one can add several 
perturbing terms and the resulting Hamiltonian are all examples of QES systems. In this 
section, we consider three such perturbing terms. 

Case I: Perturbation term. Vi = with -1 < / < 0. 

' ^ smn (X) 

The combined perturbed system is described by the Hamiltonian, 

Hj=p^ + [acosh(2a;) - M]^ + -^í^^il ^ (4.1) 

smh (x) 

We show that for integral valúes of M — Z — 1 or M + Z, this is a QES problem and one 
can obtain p QES eigenstates in case M = / + l+ p + sorM = —l+p + s respectively, 
where p = 1,2, 3, ... and s = 0, 1. 

We substitute 

^(x) =e-t^°^'^(2"V(a;), (4.2) 
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in the Schroodinger equation Hil) = Eip with H as given by Eq. (14. ip and obtain 



(j)"{x) - 2asinh(2x)0'(x) + [{E - - 

+2(M - l)a cosh(2x) - ^^\\}\ ](j){x) = O (4.3) 

sinh (x) 

On further substituting 

= [sinh(x)]°?7 , (4.4) 



we obtain 



r]"{x) + 2[a coth(a;) - a smh{2x)]r]' {x) + [E - M'^ - 

W - 2(M - l)a + 4(M - a - l)a cosh'^{x)]ri{x) = O (4.5) 



provided 



a(a-l) = /(/ + !). (4.6) 

Eq. (14.61) implies either a = / + 1 or a = — /. We first consider a = l + 1 and then it is 
easy to see that Eq. (14.51) with a = l + 1 has p QES solutions in case M = l+ l + p+ s 
where s = O or s = 1. In particular the solutions are of the form 

oo 

7] = [cosh(x)]" Y. an[cosh\x)r , (4.7) 

n=0 

where s = O in case M = l + 2p while s = 1 in case M = / + 2p + 1. Few low lying 
solutions are 

M = 1 + 2, 7] = constant , E = o? + 2a{l + 1) + 2/ + 3 , 

M = / + 3 , r] = cosh(x) , E = + 2al + 21 + b , (4.8) 

M = / + 4, E = ?/ + a2 + 3(2/ + 5) + 2(/ + 3)a, 

T] = Acoúi^{x)+ B , y = -2{l + 2 + 2a) ± 2J{1 + 2 + 2af - 4a . (4.9) 



M = / + 5, E = í/ + a2 + 8(/ + 3) + 2(/ + 4)a, 
T] = A cosh^(x) + B cosh(x) , 



y = -(Al + 9 + 8a) ± 2i/(/ + 2a)2 + 3(2/ + 3) . (4.10) 

The results for a = —l are immediately obtained from above by replacing everywhere 
/ by -/ - 1. 
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By making the substitution cosh^(x) = í in Eq. ( 14. Sp . one can show that the corre- 
sponding Bender-Dunne polynomials satisfy three term recursion relation. 

Case II: Perturbation term, Vo = — ^^rir^- 

' ^ cosh^ (x) 

The Hamiltonian of the system is thus given by 

iÍ2=p' + [acosh(2x)-M]2-i^^-tlI, (4.11) 

cosh [x) 

where / is any real number. We again show that for integral valúes of M — / — 1 or M + /, 
this is a QES problem and one can obtain p QES eigenstates in case either M = l + l+p+s 
OT if M = —l + p + s respectively, where p = 1,2, 3, ... and s = 0, 1. 

We substitute 

^(x) = e-t '^"^^'(^x) cosh''{x)T], (4.12) 
in the Schroodinger equation H2'ip = Eíp with H2 as given by Eq. fl2.14p to obtain 

r]"{x) + 2[atanh(a;) - a sinh(2a;)]V(a;) + [E - - 

+2a(2a - M + 1) + + 4(M - a - 2)a cosh^ {x)]ri{x) = O , (4.13) 

provided 

a(a-l) = /(/ + l). (4.14) 

This implies a = / + 1 or a = — /. let us consider first, a = l + 1. It is easy to see that 
Eq. f l4.13p with a = l + 1 has p qes solutions in case M = l + l+ p + s where s = O or 

s = 1. In particular the solutions are of the form 

00 

7] = [sinh(x)]" an[smh\x)]'' , (4.15) 

n=0 

where s = O in case M = l + 2p while s = 1 in case M = / + 2p + 1. Few low lying 
solutions are 

M = 1 + 2, T] = constant , E = - 2a{l + 1) + 2/ + 3 , 

M = l + 3, eía = sinh(x), E = - 2al + 21 + 5 , (4.16) 



M = / + 4, E = y + a^ + 3(21 + 5) -2a{l-l), 
1] = A cosh^(x) + B , 

y = -2{l + 2 + 2a) ± 2^(/ + 2 + 2a)2 - 4a(2/ + 3) 



(4.17) 
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M = / + 5, E = y + a^ + 8{l + 3)-2a{l-2) 
í] = A smh^(x) + B sinh(x) , 



y = -[Al + 9 + 4a) ± 2^(/ - 2a)2 + 3(2/ + 3) . (4.18) 

It is easy to convince oneself that all the above solutions are still solutions with re- 
placement of / by — / — 1 everywhere. 

The three term recursion relations for the associate polynomials is obtained in this 
case by making the substitution sinh^(x) = í in Eq. fl4.13p . 

Case III: Perturbation term, = - ^^t^. 

' ^ smn (x) cosn (x) 

We now show that it is still a QES problem even if the perturbation is the sum of the 
two perturbations considered in Case I and Case II, i.e. consider the Hamiltonian 

H,=p^ + [acosh(2x) - Mr + - ^^^^ , (4.19) 

smh [X) cosh [x) 

where — 1 < / < O so that the singularity at a; = O is not strong enough. We show that for 
positive integral valúes of either (M - / - 51 - l)/2 or (M + Z + 51 + l)/2 or (M + / - 5()/2 
or (M — l + g)/2, this is a QES problem and one can obtain p QES eigenstates in case 
either M = l + g + 2p+ 1 oí M = -l - g + 2p - 1 ov M = -l + g + 2p oí M = l - g + 2p, 
where p = 1, 2, 3, .... 

We substitute 

i¡j{x) = e-^^°'^^'^''\cosh{x)f[smh{x)]'^r], (4.20) 

in the Schroodinger equation H^ip = Eip with as given by Eq. fl4.19p we find that QES 
solutions exist only when a = g + 1, —g and /3 = / + 1, — /. On choosing a = g + 1, P = l + l, 
we obtain 

r]"{x) + 2[{l + 1) coth(x) + {g + l) tanh(x) - asinh(2x)]?7'(x) 

+ y + zcosh^ix) r]{x) = 0, (4.21) 

where 

y = E - - + {l + g + 2f + 2a{2g - M + 3) , z = Aa{M - l - g - 3) . (4.22) 

It is easy to see that Eq. fl4.2ip has p QES solutions in case M = l + g + 2p+l with 

p = 1,2,3,.... In particular the solutions are of the form 

00 

77= ^P„[cosh2(x)]", (4.23) 

n=0 
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in case M = l + g + 2p + l. Two of the low lying solutions are 

M = 1 + g + 3, r] = constant , E = - 2a{2g + 3) + 2/ + 2^ + 5 , (4.24) 

M = l + g + 5, T] = Acos\í^{x) + B , 

y = -2{l + ^ + 3 + 2a) ± 2^(/ + ^ + 3 + 2af - Aa{2g + 3) . (4.25) 

The results for the remaining three cases are immediately obtained from here by replacing 
(Z, g) with (— / — 1, —g — 1), or with (— / — 1, g) or with (/, —g — 1). 

V. CONCLUSION 

In this paper, we have shown that the known QES DSHG (and henee DSG) systems 
(both Hermitian and complex PT-invariant one) can be further enlarged by adding per- 
turbations and still it continúes to be a QES system. It will be interesting to look at other 
QES examples and obtain new QES systems by adding suitable perturbating terms. 
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